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There are at least three cosmic backgrounds of primordial gravitational waves coming 
from inflation: those produced during inflation and associated with the stretching of quantum 
T-H modes; those produced at the violent stage of preheating after inflation; and those associated 

with the self-ordering of Goldstone modes if inflation ends via a global symmetry breaking 
scenario, like in hybrid inflation. Each background has its own characteristic spectrum with 
^ specific features. We describe in detail the origin of extra peaks in the PGW background 

from preheating in the case in which inflation ends with the formation of cosmic strings due 

Qto an Abelian Higgs model after hybrid inflation. We then discuss the prospects for detecting 
each gravitational wave background, and distinguishing between them with a very sensitive 
probe, the local B-mode of the cosmic microwave background polarization. 

§1. Introduction 

Cosmological InflatioiIlI''lHl naturally generates a spectrum of density fluctua- 
tions responsible for large scale structure formation which is consistent with the 
observed CMB anisotropiesHf It also generates a spectrum of gravitational waves, 
O whose amplitude is directly related to the energy scale during inflation and which 

induces a distinct B-mode polarization pattern in the CMBllI Moreover, Inflation 
typically ends in a violent process at preheatingp'' where large density waves collide 
at relativistic speeds generating a stochastic background of GW* ' with a non-thermal 
I spectrum characterized by a prominent peak at GHz frequencies for GUT-scale mod- 

els of inflation (or at mHz-kHz for low scale models of inflation), and an amplitude 
proportional to the square of the mass scale driving/ending inflation. Such a back- 
O ground could be detected with future GW observatories like Adv-LIGO,^' LISApf 

<N BBO^ DECIGOllol etc. In case the end of inflation involves the presence of gauge 

CN fields associated with the breaking of some symmetry, like in most scenarios of hy- 

brid inflation, then the PGW spectrum presents an extra peak associated with the 
mass scale of the corresponding gauge fields 113 This would constitute a very clear 
^ signature of the physics of reheating, and its detection would open a new window 

■'-j into the early universe. 

rN Furthermore, if inflation ended with a global 0(N) phase transition, like in 

certain scenarios of hybrid inflation, then there is also a GWB due to the contin- 
uous self-ordering of the Goldstone modes at the scale of the horizonjUJ which is 
also scale- invariant on subhorizon scales JiH with an amplitude proportional to the 
quartic power of the symmetry breaking scale, that could be detectable with laser 
interferometers as well as indirectly with the B-mode polarization of the CMB ED' 
Gravitational waves produced during inflation arise exclusively due to the quasi- 
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Fig. 1. The time evolution of the GW energy density during the initial stages of preheating after 
hybrid inflation, from Ref. [6]. Note the three stages of tachyonic growth, bubble collisions and 
turbulence. 



exponential expansion of the Universe,0 and are not sourced by the inflaton fluctua- 
tions, to first order in perturbation theory. They have an approximately scale invari- 
ant and Gaussian spectrum whose amplitude is proportional to the energy density 
during inflation. GUT scale inflation has good chances to be discovered (or ruled 
out) by the next generation of CMB anisotropies probes, Planclilll COrT^P^ and 
CMBpolll3At the end of inflation, reheating typically takes place in several stages. 
There is first a rapid (explosive) conversion of energy from the inflaton condensate to 
the fields that couple to it. This epoch is known as preheating^' and occurs in most 
models of inflation. It can be particularly violent in the context of hybrid inflation, 
where the end of inflation is associated with a symmetry breaking scenario, with a 
huge range of possibilities, from GUT scale physics down to the Electroweak scale. 
Gravitational waves are copiously produced at preheating from the violent collisions 
of high density waves moving and colliding at relativistic speedspf see Fig. 1. The 
GW spectrum is highly peaked at the mass scale corresponding to the symmetry 
breaking field, which could be very different from the Hubble scale. 

In low-scale models of hybrid inflation it is possible to attain a significant GWB 
at the range of frequencies and sensitivities of LIGO or BBO. The origin of these GW 
is very different from that of inflation. Here the space-time is essentially static, but 
there are very large inhomogeneities in the symmetry breaking (Higgs) field due to 
the random spinodal growth during preheating. Although the transition is not first 
order, "bubbles" form due to the oscillations of the Higgs field around its minimum. 
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The subsequent collisions of the quasi-bubble walls produce a rapid growth of the 
GW amplitude due to large field gradients, which source the anisotropic stress- 
tensorHf The relevant degrees of freedom are those of the Higgs field, for which 
there are exact analytical solutions in the spinodal growth stage, which later can be 
input into lattice numerical simulations in order to follow the highly non-linear and 
out-of-equilibirum stage of bubble collisions and turbulence. However, the process 
of GW production at preheating lasts only a short period of time around symmetry 
breaking. Soon the amplitude of GW saturates during the turbulent stage and then 
can be directly extrapolated to the present with the usual cosmic redshift scaling. 
Such a GWB spectrum from preheating would have a characteristic bump, worth 
searching for with GW observatories based on laser interferometry, although the 
scales would be too small for leaving any indirect signature in the CMB polarization 
anisotropics. Moreover, the mechanism generating GW at preheating is also active 
in models where the SB scenario is a local one, with gauge fields present in the 
plasma, and possibly related to the production of magnetic field fiux tubes JHI' In 
such a case, one could try to correlate the GWB amplitude and the magnitude and 
correlation length of the primordial magnetic field seed. 

§2. Gravitational waves from Abelian Higgs Cosmic Strings 

Gravitational waves are a robust prediction of general relativity. There is indirect 
evidence of their existence from inspiraling binary pulsars, although no single direct 
detection has been claimed. A stochastic background of GW may soon be discovered, 
either directly with laser interferometer antennas or indirectly through the pattern 
of polarization anisotropics they induce in the cosmic microwave background. Such 
a detection would open a completely new and unexplored window into the Early 
Universe, possibly as rich as that which has been recently revealed in the CMB. 
There are many sources of GW that can generate a stochastic backgrounds and thus 
it is necessary to characterize those backgrounds with as much detail as possible. 
Apart from known astrophysical point-like sources beyond the confusion limit (where 
we cannot resolve them) , there are also predictions for GWB from cosmic defects and 
hypothetical strongly first order phase transitions in the Early Universe. Moreover, 
Cosmological Inflation makes a robust prediction of a stochastic GWB produced 
during the quasi-exponential expansion of the Universe, with very specific spectral 
signatures: a Gaussian, almost scale-invariant spectrum with an amplitude directly 
related to the energy scale of infiation. If the scale is high enough (close to the 
Grand Unification Theories' scale) then these GW will also leave an imprint in the 
(curl) polarization anisotropics of the CMB. Unfortunately, this GWB is still too 
weak to be discovered with the near-future GW interferometric antennas, although 
it does cover a sufficiently broad frequency range to be detectable by future GW 
Observatories (GWO) like BBO or DECIGO. 

Furthermore, a robust prediction of inffation is that it must have ended, convert- 
ing the vacuum energy responsible for the tremendous expansion into the matter and 
radiation we observe today. Such a process, known as reheating, is typically very vio- 
lent and very inhomogeneous, with large density waves moving at relativistic speeds 
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and colliding among each other, thus converting a large fraction of their gradient 
energy into gravitational waves. In some cases, the conversion is so sudden and vi- 
olent that a significant fraction of the total energy that goes into radiation ends in 
a stochastic background of GW, which could be detected in the future. The energy 
spectrum of such a GWB is very non-thermal and far from scale invariant, but ac- 
tually peaked at a frequency which is related to the typical mass scale responsible 
for the end of inflation (either the mass of the inflaton or that of the field that trig- 
gers the end of inflation, like in hybrid models), which could be orders of magnitude 
smaller than the Hubble scale at the end of inflation. However, if the energy scale of 
inflation is large (of order the GUT scale) then this stochastic GWB will be peaked at 
GHz frequencies, far from the present sensitivity of GW interferometers. Nowadays, 
our only chance of detecting the GWB from reheating is to consider the low-scale 
models of inflation - like hybrid models - with the appropriate parameters to convert 
a large fraction of the initial vacuum energy into GW. The analyses done so far have 
considered only scalar fields whose gradient energies source the anisotropic stresses 
needed for GW production. However, vector fields (gauge or not gauge) are expected 
to be an even better source of GW, due to their anisotropic curl components, so that 
preheating scenarios with gauge fields may have a larger contribution to the GWB 
than scalar models. In fact, previous studies of gauge fields at preheating, in the 
context of Electroweak Baryogenesis and in the generation of Primordial Magnetic 
Fields, have identified the formation of long- wave semiclassical gauge field config- 
urations like sphalerons and helical magnetic fiux tubes which evolve with time very 
anisotropically and could contribute significantly to the production of GW. 

In a recent paper jI3 we developed a formalism to calculate the production of GW 
by coupled systems of scalar and gauge fields on the lattice. The numerical method 
that we have constructed can be applied to different sources of GW where out-of- 
equilibrium gauge fields play an important role, such as thermal phase transitions, 
cosmological networks of local defects or non-perturbative decays of scalar conden- 
sates into gauge fields. We have studied in detail the dynamics and the production of 
GW during preheating after hybrid inflation, in the context of abelian-Higgs models 
that go through dynamical symmetry breaking triggered by the expectation value of 
the inflaton fleld. As the inflaton is driven slowly (as opposed to a quench) below the 
critical value, the mass squared of the abelian Higgs fleld becomes negative and drives 
the spinodal growth of long-wave modes of the Higgs. Since the Higgs is charged, 
its rapid growth induces a corresponding growth of gauge field configurations. At 
the end of inflation there are no temperature fluctuations that can induce over-the- 
barrier transitions. However, long- wave quantum fluctuations become semi-classical 
and act as a stochastic force that allow transitions over the false vacuum and thus 
induce (locally) the generation of a topological winding number of the Higgs fleld. 
After symmetry breaking, there is not enough energy to unwind the Higgs phase, 
leaving behind a Nielsen-Olesen string. Such cosmic string conflgurations can be 
seen explicitly in our spatial distributions of both Higgs and gauge fleldsED They 
play a crucial role in the production of GW at preheating and we observe that the 
spatial distribution of GW is indeed concentrated around the strings. Those strings 
will eventually decay (we see that they become wider and disperse away their en- 
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ergy density in the form of small-scale structures of the fields, although the winding 
phase around the core of the string remains, since it is topologically stable), which 
eventually shuts-off the GW production. 
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Fig. 2. The fraction of energy density in GW as a function of the typical momenta at preheating 
after inflation in the abelian-Higgs model for different ratios of the gauge field to Higgs mass. 
Note the existence of three peaks in the most extreme cases: an infrared peak due to straight 
segments of strings, and two more ultraviolet peaks, related to the Higgs and gauge field masses. 

The complicated dynamics occurring at preheating in this abelian Higgs-inflaton 
model has been studied using both power spectra analyses, as well as the fields' dis- 
tributions in configuration space, together with histograms of the fields' values, as a 
function of time, in order to correlate the different features observed and their evo- 
lution. The picture that arises is the following. At the end of inflation the tachyonic 
growth of the Higgs Gaussian random field creates an inhomogeneous distribution 
of fields characterized by "bubbles" of Higgs energy density that expand and collide. 
The gauge field concentrates at the valleys between the bubbles, where the Higgs has 
low values, forming long flux tubes of magnetic energy density. The dynamics of the 
bubles when they expand and collide leads to regions in space where the Higgs field 
reaches the false vacuum and there are over-the-barrier transitions, with topological 
windings associated with them. These Nielsen- Olesen vortices are connected with 
each other in a cosmic string which runs along the core of the magnetic flux tubes. 
There are strings that encompass the whole simulation box and even beyond, thanks 
to periodic boundary conditions. We observe this process both in configuration space 
and with the histograms of Higgs vevs. 

We have followed the dynamics of the strings during and after the symmetry 



6 



Juan Garcia-Bellido 



breaking, although still on time scales shorter than the Hubble time and on length 
scales smaller than the Hubble radius. Once the strings are formed, they evolve by 
increasing their size and shedding away layers of magnetic energy density. At the 
cores of the strings there always remain a thin magnetic flux line but the energy 
seems to pour away from the strings in the form of waves concentric with the string. 
Nevertheless, we observe (in a transverse plane to the string) that at the core of the 
string there remains a conserved winding number of the Higgs. We have followed 
this winding number up to long times and we confirm that it is still there, in spite 
of the fact that the magnetic flux tube is so dilute that we cannot see it coherently: 
it seems to have "evaporated" . What remains is a diffuse background of small-scale 
structures of the Higgs and gauge fields permeating the whole box, together with 
the remnants of the strings. 

The formation, evolution and fragmentation of the strings are accompanied by 
a significant production of gravitational waves which inherit specific features from 
the string dynamics. In position space, we observe how the distribution of GW 
follows very closely the evolution of the strings, being first concentrated around the 
straight segments of strings, then fattening as the strings become wider and finally 
being dispersed over the lattice as the strings emit small-scale structures of the fields. 
In Fourier space, this dynamics is encoded into the successive appearance of very 
distinct peaks in the GW spectra. The position of each peak is directly related to 
the physical scales in the problem: the Higgs mass, which governs the width and 
interactions of Higgs field's strings, the gauge field mass, which governs the width 
and interactions of gauge field's strings, and the typical momentum amplified by 
tachyonic preheating, which determines the characteristic size of the bubbles when 
they collide and the correlation length of the straight segments of strings. The former 
two determine the peaks in the high momentum (UV) range of the spectrum, while 
the latter corresponds to the long-wave (IR) peak. The IR peak appears first, when 
the bubbles collide and the strings are formed, while the UV peaks are formed later 
on, when the strings evolve and decay into small-scale structures of the Higgs and 
gauge fields. When the different scales are close to each other, the different peaks 
are superimposed and the amplitude in GW increases. When the gauge coupling 
constant is significantly smaller than the Higgs' self-coupling, the results reduce to 
the GW spectra produced without gauge field, characterized by a single peak. 

We have calculated the GW spectra produced in this abelian Higgs model of 
preheating after hybrid inflation with state-of-the-art simulations, although still lim- 
ited in spatial resolution and box sizes. In order to probe reliably the different scales 
in the problem in each simulation, we developed a lattice calculation of GW produc- 
tion with gauge fields that is accurate up to second order in the lattice spacing. The 
present-day frequency and amplitude of these GW are very sensitive to the model 
parameters and the frequency of the different peaks may differ by many orders of 
magnitude, as illustrated in Fig. [2] As in the same model with only scalar fields, very 
small coupling constants are still neccessary for these GW to fall into a frequency 
range that is accessible by interferometric experiments. Whether this is natural or 
not depends on the underlying theory for inflation and particle-physics models of hy- 
brid inflation with such small coupling constants have indeed been already proposed 
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in the literature, sed^} and references therein. We also observed that the frequency of 
the IR peak in the GW spectrum can be smaller than the peak frequency produced 
in the same model with only scalar fields, so the gauge field may enlarge the regions 
of the parameter space that may lead to an observable signal. More generally, there 
are many other models of inflation and preheating where gauge fields may play an 
important role and wich may lead to GW that could be observed in the future. 




10-6 0.001 1 1000 10^ 10'' 

f(Hz) 

Fig. 3. The predicted stochastic background of GW from preheating in the abehan-Higgs model for 
two different sets of parameters, together with the expected sensitivity from future GWO hke 
Advanced LIGO/VIRGO, LISA, ET, BBO and DECIGO. Note that in order to reach GWO 
sensitivity we had to extrapolate the position of the IR peaks and make an educated guess for 
the shape of the spectra between the peaks (dashed hues) . Also plotted are the expected GWB 
from a global phase transition and that from an inflationary model with tensor to scalar ratio 
r = 0.1. 

After preheating the system enters into a turbulent regime where, at least in the 
abelian Higgs model that we considered, gravitational waves are no longer produced 
and the GW energy density saturates. We expect this result to be rather generic 
in abelian scalar gauge theories, because in that case the gauge fields that are pro- 
duced at preheating acquire a mass, either directly through the Higgs mechanism or 
due to their interactions with scalar fields' fiuctuations. It would be interesting to 
study models with other gauge groups, like SU{2) x U{1), where gauge fields remain 
effectively massless after symmetry breaking. This may happen for instance during 
preheating after hybrid infiation close to the electroweak scale (possibly a secondary 
stage of infiation, not necessarily related to the CMB anisotropics, only responsible 
for reheating the universe), where the photon spectra may exhibit inverse cascade 
during the turbulent evolution towards thermal equilibrium. This could signifi- 
cantly lower the typical frequency of the resulting GW today and relax the conditions 
on the parameters for these GW to be observable. The details of the GW spectra 
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produced from preheating should also be rather sensitive to the particular gauge 
group under consideration because this determines the nature of the defects that 
can be formed. The defects do not have to be stable since these GW are produced 
when they are being formed. 

An intriguing possibility is the following. Given that preheating is so extremely 
inhomogeneous, and since these inhomogeneities get imprinted in the gravitational 
wave background, which immediately decouples from the plasma, one may envision 
a stage of technological development in the not so far future in which GWO with 
sufficient angular resolution may resolve the structures that gave rise to the GWB 
right at the moment when the Universe reheated. In the usual preheating scenario at 
high energy scales with only scalar fields, the physical structures will have today a size 
that is completely undetectable when projected over the sky, and thus the GWB will 
look essentially homogeneous from Earth. However, if gauge fields with long string- 
like configurations (of horizon size and possibly even with superhorizon correlations) 
were behind the generation of the GWB, then one could expect to see inhomogeneities 
in the angular distribution of those gravitational waves. In particular, an array of 
GWO could detect the string-like anisotropics in the GWB across the sky. At the 
moment, the angular resolution of LIGO is not better than a degree projected in 
the sky. However, in the future one could resolve much finer structures in the GWB 
thanks to a dense network of ground-based laser interferometersEB' Thinking ahead 
of our times, it may not be unrealistic to imagine that in the not so far future the 
GWB will be mastered with sufficient detail to resolve the anisotropics in this elusive 
background and thus recover vital information about the physics responsible for the 
violent conversion of energy from inflation to a radiation and matter dominated epoch 
(the Big Bang of the Old Theory). No other probe can give us so much information, 
since GW decouple immediately upon production and thus retain the spatial and 
energy distributions of the sources that produced them. We can compare with the 
Cosmic Microwave Background, which gives us detailed information about the epoch 
of photon decoupling thanks to the exquisite measurements of the angular correlation 
of both temperature and polarization anisotropics. The CMB provides a snapshot 
of how the Universe was like 380,000 years after the Big Bang. On the other hand, 
the GWB would open a window into the physics of the Big Bang itself, allowing 
us to infer from its detailed features whether it was as violent and inhomogeneous 
as we predict, to determine what kind of fields were present and whether the rich 
phenomenology that we associate with preheating (topological defects, baryogenesis 
and/or leptogenesis, primordial magnetic seed creation, non-thermal production of 
dark matter, etc.) was actually realized in nature. 

§3. Gravitational waves from global defects formed after inflation 

In the case that inflation ends with a global or local symmetry breaking mech- 
anism, then there generically appear cosmic defects associated with the topology of 
the vacuum manifold. For instance, global cosmic strings are copiously produced 
during preheating if the Higgs field is a complex scalar with a U(l) global symme- 
trylH In principle, all kinds of topological and non-topological defects could form 
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Fig. 4. The local B-polarization power spectra for tensor perturbations from inflation, cosmic 
strings, textures and the large- TV limit of thenon-linear sigma-model. All spectra are normalized 
such that they make up 10% of the observed temperature anisotropy at ^ = 10. Note that the 
lensed EE modes (red dashed line) contribute as colored noise to the local-BB angular power 
spectrum. 



at the end of inflation and during preheating. Such defects will have contributions 
to all three different metric perturbations: scalar, vector and tensor, with similar 
amplitudesJl^} In a recent workjUl' we analyzed the production of GW at preheating 
for a model with 0(N) symmetry. The dynamics at subhorizon scales was identical 
to that of the usual tachyonic preheating. However, in this model even though the 
Higgs potential fixed the radial component to its vev, there remained the free (mass- 
less) Goldstone modes to orient themselves in an uncorrelated way on scales larger 
than the causal horizon. In the subsequent evolution of the Universe, as the horizon 
grows, spatial gradients at the horizon will tend to reorder these Goldstone modes in 
the field direction of the subhorizon domain. This self-ordering of the fields induces 
an anisotropic stress-tensor which sources GW production. In the limit of large N 
components, it is possible to compute exactly the scaling solutions, and thus the 
amplitude and shape of the GWB spectrum. It turns out that the GWB has a scale- 
invariant spectrum on subhorizon scaleflf and a infrared tail on large scalesjlD 
which can be used to distinguish between inflation and these non-topological defects. 

Apart from the IR tail, the main difference between inflationary and global defect 
contributions to the CMB anisotropies arises from the fact that defects generically 
contribute with all modes: scalar, vector and tensor modes, with similar amplitudes, 
while inflationary tensor modes could be negligible if the scale of inflation is well 
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Fig. 5. The local B-polarization angular correlation functions for 6 < 1° for inflation, cosmic 
strings, textures and the large- A'^ limit of NL sigma model. 

below the GUT scale. Since (curl) B-modes of the polarization anisotropics only 
get contributions from the vector and tensor modes, the detection of the B-mode 
from inflation may be challenging, and dedicated experiments like Planck, COrE 
and CMBpol have been designed to look for them. On the other hand, defects' 
contribution to the temperature anisotropics have a characteristic smooth hump in 
the angular power spectrum, which allows one to bound their amplitude (and thus 
the scale of symmetry breaking) below 10^^ GeV.'^ However, the contribution to 
the B-mode coming from defects have both tensor and vector components, and the 
latter can be up to ten times larger than the former, and actually peaks at a scale 
somewhat below the horizon at last scattering (in harmonic space the corresponding 
multipole is £ ~ 1000). 

In a recent papeilHl we analyzed the possibility of disentangling the different 
contributions to the B-mode polarization coming from defects versus that from in- 
flation. The main difficulty, for both defects and tensor modes from inflation, is 
that the B-mode power spectrum is "contaminated" by the effect of lensing from 
the intervening matter distribution on the dominant E-mode contribution on similar 
angular scales. Using the temperature power spectrum to determine the underlying 
matter perturbation from evolved large scale structures responsible for CMB lens- 
ing, it is possible to ejigineer an iterative scheme to clean the primordial B-modes 
from lensed E-modesH?! This procedure leaves a significantly smaller polarization 
noise background /\p,eff which allows one to detect the GW background at high 
confidence level (3-ct) if the scale of inflation or that of symmetry breaking is high 
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Fig. 6. The signal-to-noise ratio as a function of the normahzed polarization sensitivity, for inflation, 
cosmic strings, textures and the large- A'^ limit of the non-linear sigma-model. 

enough. What we reahzed is that the usual E- and B-modes used for computing the 
angular power spectra are complicated non-local functions of the Stokes parameters, 
involving both partial differentiation and inverse laplacian integration. Such a non- 
local function requires knowledge of the global polarization on scales as large as the 
horizon, where the B-mode angular correlation function is negligible and thus prone 
to large systematic errors. In contrast, the so-called "local" E- and B-modeaU'lMl' 
can be constructed directly from the Stokes parameters and do not involve any non- 
local inversion. A direct consequence in this change of variables is that the angular 
power spectrum of local B-modes has a extra factor ni = (£-|-2)!/(-£ — 2)! ~ which 
boosts the high-£ peak in the defects' power spectra. 

When compared with the angular correlation function of inflation, it gives a sig- 
nificant advantage to the defects' prospects for detection in future CMB experiments, 
see Fig. [7} In Table 1 we give the expected bounds at 3-a on the tensor amplitude 
for inflation and various defect models, while in Table 2 we give the corresponding 
bounds on the scale of inflation, and on the symmetry breaking scale for defects, 
according to 

Minf = 1.63 X 10^6 GeV (^^) , v = 1.02 x 10^^ GeV 
This gives the fundamental connection between the CMB observations and the par- 
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Fig. 7. The relative signal-to-noise ratio for defect models versus inflation for local (continuous 
lines) and non-local (dashed lines) B-modes. 



tide physics models at ultra high energy, much beyond any experiment will ever be 
able to reach. 

Table I. The limiting amplitude for inflation (r=T/S) and various defects (e — Gv^), at 3-a in the 
range 6 G [0, for Planck {Ap,c« ~ 11.2 /xK-arcmin), COrE-like exp. {Ap^cs ~ 0.7 ^K-arcmin) 
and a dedicated CMB experiment with Ap^^a ~ 0.01 /iK-arcmin. We have assumed /sky = 0.7 
for all CMB experiments. 
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